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ABSTRACT 

The  Griffith -Irwin  (or  linear -elastic)  fracture  concept  was  used  to  investigate 
the  feasibility  of  the  application  of  fracture  mechanics  to  orthotropic  plates.  The 
problem  of  an  infinite  plate  containing  a  single  crack  orientated  in  the  direction  of  one 
of  the  planes  of  elastic  symmetry  of  the  plate  subjected  to  arbitrary  plane  loading  was 
examined.  Two  analytical  crack  tip  stress  analyses  for  this  case  were  presented, 
and  unlike  the  isotropic  case,  the  stress  distributions  were  observed  to  be  dependent 
on  the  material  constants.  However,  as  in  the  isotropic  case  the  crack  tip  stress 
singularity  was  observed  to  be  r  '  and  stress  intensity  factors  k^  =  crV a  and  k^  =xV~a 
for  loading  symmetric  and  skew -symmetrical  to  the  crack  respectively  could  be  used. 
Further,  the  existence  of  a  functional  relation  between  k^  and  k 2  was  proposed. 

Results  of  the  experimental  investigation  using  balsa  wood  plates  indicated  that 
the  critical  stress  intensity  factors  k  and  k2C  were  constants  for  tension  and  pure 
shear.  Moreover,  under  combined  tension  and  shear,  klc  and  k2c  were  found  to  be 

k  1  k  2 

related  by  the  function  (-rr— )  +  (t - — )  =  1  which  satisfied  the  physical  considerations. 

Klc  K2c 
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I.,  PRELIMINARY  REMARKS 

In  many  engineering  applications,  flaws  and  minute  cracks  are  observed  to  extend, 
leading  to  fracture  at  nominal  stress  levels  below  that  of  the  yield  strength  of  the  struc¬ 
tures.  The  study  of  the  mechanics  of  fracture  and  attempts  to  predict  the  combination 
of  crack  length  and  stress  which  would  lead  to  failure  by  spontaneous  crack  extension 
have  been  undertaken  by  many  investigators.  Fracture  mechanics  thus  developed  in  the 
past  20  years  has  achieved  remarkable  success  in  applications  to  isotropic  materials. 
Since  many  engineering  materials  are  anisotropic,  it  is  logical  to  analyze  the  parameters 
involved  and  examine  whether  fracture  mechanics  can  be  applied  to  anisotropic  materials 
as  well. 

A  special  but  important  group  of  anisotropic  materials,  the  orthotropic  material, 
was  chosen  for  analysis.  This  investigation  is  applicable  to  single  as  well  as  composite 
materials  of  two  or  more  components  which  are  grossly  homogeneous,  but  which  possess 
orthotropic  properties.  Some  important  examples  are  wood,  oriented  glass  fiber  rein¬ 
forced  plastics  and  high  strength  wire  reinforced  metals  and  plastics. 

The  investigation  consists  of  two  parts.  In  the  first  part,  a  theoretical  analysis 
of  the  stress  distribution  in  the  vicinity  of  the  tip  of  a  crack  in  an  orthotropic  plate  is 
presented.  In  the  second  part,  an  experimental  program  is  designed  to  verify  the 
analytical  findings.  Thi6  study  of  orthotropic  materials  was  undertaken  in  an  attempt 
to  shed  light  upon  the  application  of  fracture  mechanics  to  general  composite  materials 
made  of  metals  as  well  as  non-metals. 
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II.  ANALYSIS 

A.  Introduction  to  Analysis 

The  quasi-static  stability  criterion  of  a  crack  in  a  linear  isotropic  material 
was  formulated  by  Griffith  [l]  *  and  later  extended  by  Irwin  [2j  .  In  its  simplest  form 
the  linear-elastic  (or  Griffith-Irwin)  fracture  theory  states  that  the  critical  energy  dissi¬ 
pation  rate  per  unit  area  associated  with  an  unstable  crack  extension  is  a  constant  for  a 
given  material.  In  the  formulation  of  this  theory,  the  following  assumptions  are  made: 
1)  the  material  is  macroscopically  isotropic  and  homogeneous,  2)  the  material  is  ideally 
elastic,  3)  infinitesimal  linear  elasticity  is  applicable  for  stress  analysis  in  the  vicinity 
of  the  tip  of  the  crack.  Irwin,  upon  defining  three  kinematically  possible  modes  of 
crack  extension,  was  able  to  compute  the  strain  energy  release  rates,  ^  2  and  y 

associated  with  each  mode. 

Further,  in  examining  the  elastic  stress  analysis,  it  was  observed  that  the  stress 
intensity  factors  k  ,k„  and  kq  for  each  of  the  three  modes  are  directly  related  to  the  re- 
spective  strain  energy  release  rates  €j  ^  ,  "Cj  2  and  Cj y  Consequently,  the  magnitude 
of  the  critical  stress  intensity  factors  was  also  associated  with  unstable  crack  propa¬ 
gation.  Based  on  the  linear -elastic  fracture  theory,  many  problems  in  brittle  fracture 
of  isotropic  materials  have  been  described  with  a  useful  degree  of  engineering  accuracy 
when  similitude  was  maintained  between  the  model  and  the  actual  crack  extension  prob¬ 
lems. 

Obviously,  the  Griffith-Irwin  fracture  concept  is  not  limited  to  isotropic  mate¬ 
rials.  In  order  to  gain  insight  to  the  feasibility  of  application  of  fracture  mechanics  to 
anisotropic  materials,  the  simplest  case  of  anisotropy  -  the  orthotropic  material-is  con¬ 
sidered. 

A  material  is  said  to  be  orthotropic  when  it  possesses  three  orthogonal  planes  of 
elastic  symmetry.  The  majority  of  the  orthotropic  engineering  materials,  such  as  ply¬ 
wood  and  some  fiberglass  reinforced  plastics,  possess  such  properties.  However,  if 
the  analysis  is  limited  to  materials  which  possess  only  two  orthogonal  planes  of  elastic 
symmetry,  the  complexity  of  the  mathematical  formulation  of  their  constitutive  equations 
is  considerably  reduced.  Limiting  the  analysis  to  such  a  material  corresponds  to  the 
study  of  only  one  layer  of  the  plywood,  or  only  one  layer  of  the  reinforced  plastic  with 
glass  fibers  orientated  all  in  one  direction.  Nevertheless,  in  order  that  the  fracture 
mechanics  aspects  of  the  problem  shall  not  be  subdued  by  the  algebric  complexity  of  the 


*  Numbers  in  brackets  designate  References  listed. 
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more  general  problem,  the  use  of  a  simplified  model  has  to  be  tolerated. 

An  orthotropic  plate  which  possesses  two  elastic  symmetries  containing  a  single 
crack  subjected  to  arbitrary  plane  tractions  with  bending  absent,  as  shown  in  Fig.  lb, . 
will  be  considered.  The  imposed  conditions  are: 

1)  the  material  is  linearly  elastic 

2)  the  plate  is  thin  and  extends  infinitely  in  all  directions  from  a  crack 

3)  the  crack  is  straight  and  "sharp"  i.e. ,  a  straight  line  of  discontinuity 
in  the  plate 

4)  the  crack  is  parallel  to  one  of  the  principle  planes  of  elastic  symmetry 

5)  the  influence  of  environment  is  neglected 

Although  the  above  conditions  are  necessary  for  the  mathematical  analysis  which  follows 
they  do  not  significantly  restrict  its  applications.  If  a  crack  is  small  in  the  actual 
structure  compared  to  its  surrounding  dimensions,  a  sufficient  similitude  exists  between 
the  mathematical  model  and  the  actual  crack. 

B.  Analysis 

When  a  body  containing  a  crack  is  subjected  to  external  tractions,  the  sole 
influence  of  the  external  conditions  on  the  crack  (the  influence  of  environment  is  neg¬ 
lected)  is  the  state  of  stress  around  the  crack.  If  the  crack  is  to  extend,  the  state  of 
stress  or  other  parameters  associated  with  the  state  of  stress  at  the  tips  of  the  crack 
must  reach  a  critical  value  for  the  given  material  of  which  the  body  is  composed.  Then, 
to  analyze  the  crack  stability  in  an  orthotropic  plate,  an  examination  of  the  state  of 
stress  in  the  vicinity  of  the  tips  of  the  crack  is  in  order. 

For  the  plane  problem  of  a  linear  elastically  orthotropic  material,  the  general¬ 
ized  Hooke's  law  may  be  written: 

ex  ail  CTx  +  al2  CTy  +  al6  Txy  ’ 

ey"a120x+a22c,y  +  a26  V  (1) 


2e  =  a1,;cr  +ar,,<j  +  a, ,  t  . 
xy  16  x  26  y  66  xy* 


The  compatibility  equation  for  plane  deformation  is 


8x8y 


(2) 
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In  order  to  simplify  the  constitutive  equations,  the  x  and  y  axes  can  be  oriented  in  such 
a  way  that  they  coincide  with  the  principle  planes  of  elastic  symmetry.  For  this  orien 
tation  of  axes, 


al6  a26  0  • 


(3) 


Substitution  of  Eqs.  (2)  and  (3)  into  Eq.  (1)  gives 

X)  +  a 


a2o 


d2u 


d2  a 


‘11 


9y 


T  +  a12  <- 


+ 


9x 


9y 


22 


2  2 
9  cr  9  T 


9x 


9x9y 


(4) 


When  body  forces  are  neglected  the  equations  of  equilibrium  and  compatibility  are 
satisfied  by  Airy's  stress  function  U,  and  the  stresses  can  be  expressed  as 


o  = 
x 


2  2 
9ZU  9  U  . 

— *—  >  °  and  t 

y  *y 


9y 


2 

9  U 
9x9y 


(5) 


Also,  the  material  constants  a. ^  can  be  expressed  in  terms  of  the  usual  engineering 
notations,  i.e. 


all  E 


22 


,  Vx  V 

=  "E-  ’  ai2  =  "  7E^and 


a66  G 


(6) 


xy 


Using  the  relationships  given  by  Eqs.  (5)  and  (6),  Eq.  (4)  becomes 


94U 

^7 


+  2  <- 


2G 


xy 


v  ) 

y 


94U 

~  2  2" 
9x  9y 


E 

+  (-g2-) 


4 

9  U 

97" 


=  0. 


(7) 


The  problem  therefore  reduces  to  the  solution  of  the  above  4th  order  homogeneous 
partial  differential  equation  whose  auxiliary  equation  is 


4 

s 


+  2( 


2G 


xy 


"  v  ) 

y 


+  <*i>  =o. 


(8) 


It  can  be  shown  that  Eq.  (8)  has  no  real  roots  [3j  .  Further  since  the  coefficients  are 
real  constants  (the  elastic  constants  are  real), the  roots  must  be  in  conjugate  pairs,  i.e. 
they  are  s^  s2,  5^  s2<  [4] 

If  s1  =  (a  +  i/3  )  , 

then  s2  =  (  -a  +  i/3)  t  (9) 

and  s  =  (a  -  i/3)  , 
s2  =  -  (ct  +  i£). 


5. 


Substitution  of  Eqs.  (9)  into  (8)  gives 

E, 

:1G 


2  E  1/2  E 
2$  =(-h^)  +  (-.  X 


xy 


vx>> 


,  E  1/2  E 
2«  ^  -<*?-  -vx). 

y  xy 

7 


Then,  the  general  solution  of  Eq,  (p)  is 


(10) 


u  =  Zx( £x)  +  Z2  (ty  +  Z3  (f2)  +  Z4  (t2)  (11) 

where 

£,  =  x+  s  y,  t1  =  x  +  s  y,-  e2=x+s2y  and  t2  =  x  +  s2y. 

In  order  to  obtain  the  function  U  which  satisfies  the  boundary  conditions  for  the  plate  con¬ 
taining  a  crack,  two  solutions  will  be  outlined.  First,  a  technique  employing  a  Wester  - 
gaard  type  of  stress  function  will  be  discussed.  Following  this  a  brief  outline  of  a  stress 
function  type  of  procedure  developed  primarily  by  S.  G.  Lekhnitskii*  will  be  given. 


1.  Solution  by  Westergaard  Type  of  Stress  Function: 

From  physical  consideration,  it  is  known  that  U  must  be  real.  Therefore 
Eq.  (11)  can  be  written  in  the  following  form;** 


U  =  A  Re  [za(£1)]  +  B  Im  [ZgU^]  +  C  Re  [zc(t2)j  +  D  Im  [z D(t2)J  (12) 

where 

=  x  +  s  y  =  x  +  (a  +  i0)  y  =  x  +  iy  , 

(13) 

f2  =  x  +  s2y  =  x  +  (  +  i^)  y  =  x2  +  iy2 . 

The  arbitrary  external  tractions  are  treated  in  two  separated  parts.  1)  Loading  sym¬ 
metric  with  respect  to  the  crack  and,  2)  loading  skew -symmetric  with  respect  to  the 
crack.  The  following  notations  are  adopted: 

The  complex  functions  Z.in  Eq.  (11)  for  symmetric  loading  are 

Z1  =  Zll’  Z1^2^  =  Z12 ' 


*  A  third  method  of  solution  employing  integral  equations  by  D.  D.  Ang  and  M.  L. 
Williams  is  also  available  [5]  . 

*’  This  can  be  seen  by  expanding  the  arbitrary  functions  Z.  from  Eq.  (11)  in  power 
series.  1 


6. 


And  for  the  skew -symmetric  loading,  they  are 


Z2  Z21’  Z2  ^2^  =  Z22‘ 


In  general 


Z..  =  Z.  (£.). 

ij  i  r 


The  first  subscript  refers  to  the  nature  of  the  traction,  "1”  for  symmetric,  "2"  for 
skew -symmetric.  The  second  subscript  "1"  or  "2’'  refers  to  the  complex  variable 
and  respectively. 

Westergaard  [6]  observed  that  the  stress  around  a  mathematical  crack  on  the 
axis  y  =  0  can  be  solved  by  a  single  stress  function  in  isotropic  materials.  Irwin  [7] 
later  pointed  out  that  in  specially  orthotropic  material  (a^  =  =  0)  the  counterparts  of 

a  Westergaard  type  of  stress  functions  are  as  follows*: 

a)  Symmetric  loading 

Ux  =  -y  Re  [z  1(S  1)  +  Zjttj)]  +  ^  Im  [Z^)  -  Z1(t2)|  ,  and  (14) 

b)  Skew  symmetric  loading 


i-  1,2 
J-  1.2 


where  ** 

dZ..  _  dZ..  dZ.. 

— t5=L«Z..,  =  Z..  and  -#- =  Z\. 

d£.  ij’  dtj  ij  d£.  ij 

and 

Z..  =  Z.(£.) 
ij  1'  J7 

From  the  Cauchy -Riemann  relationships  the  following  relations  can  be  observed: 

9ReZ..  9ImZ.. 

ReZ .. - * — a - 

^  8xi  9yj  i  =  1, 2 


(16) 


9ImZ..  9Re^.. 

Im  z..  =— — - 1J 


j  =  1,2 


(17) 


ij  9x. 


J 


9yi 


*  Irwin  was  instrumental  in  adopting  Westergaard-type  of  stress  functions  to  ortho - 
tropic  materials  [7]  .  Hence  the  notation  used  by  Westergaard  and  Irwin  was 
adopted  in  this  analysis. 

**  Note,  the  summation  convention  is  not  used  here. 
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Upon  noting  that 


3x 

3x, 


=  1, 


4JL  =  o 

dxl  > 


*L  = 


3x 


a 


"5y"i  T  '  3  ’ 


9>L=  1  il.  =0  =  — 

3v  u>  Sir  ft  ’ 


2  -2  8y2 


3x 


3x  _  ja_ 

3y,  3 


and  from  Eq.  (13)  and  using  the  rules  of  differentiation  to  operate  on  Eq.  (17) »  the 
following  relations  are  obtained; 


and 


3  ReZ. 


il 


3x 

3  Im  Z. 


=  ReZ.. 


il 


Hx 


=  Im  1 


3  ReZ. 


il 


3y 

3  Im  Z 


=a  Re  Zu  -|3  Im  Z.^ 


il 


3y 


3  ReZ. 


=a  Im  Zn  +  0  Re  2{V 


i2 


3x 

3  Im  Z 


=  Re  Z 


i2 


3x 

3  ReZ 


i2 


=  Im  Zi2 ' 


il 


3y 

3  Im  Z-2 

9y 


=  -a  Re  Z.2  -0  Im  Z.^ 
=  -a  Im  Zi2  +P  Re  Z.2> 


i  =  1,2 


(18) 


i  =  1,2 


(19) 


With  these  equations,  the  stress  distribution  for  the  case  where  loading  is  symmetric 
to  the  crack  and  the  case  where  loading  is  skew -symmetric  to  the  crack  can  be  computed, 
a.  Case  l;  Symmetric  Loading 

Returning  to  Eq.  (14)  and  adopting  the  notation  which  was  described 
following  Eq.(l3),  the  Westergaard -Irwin  stress  function  can  be  written  as 

Ul  =  -lRe  Zn  +  f12  Zll‘Z12  • 


(20) 
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Using  the  results  of  Eq.  (18)  and  (19)  for  i  =  1,  substitution  of  Eq.  (20)  into  Eq.  (5) 
gives 

4^  R'<Zn+Z12>-£  Im  <ZU  '  Z12>  • 

<’y'TRe<ZU  +  Z12>+  4  Im  (Z11-Z12>' 


(21) 


2,a2 

T  =  -  “±JL- 

XV  lot 


Re  (Zu  -  Z12) 


For  an  infinite  plane  containing  a  single  crack  extending  from  x  =  a  to  x  =  -  a,Wester- 
gaard  [6]  showed  that  the  function  Z  (£)  =  satisfied  the  boundary  conditions 

and  produced  a  branch  cut  on  -  a<x<  +  a.  Since  no  consideration  of  the  material 
properties  of  the  plate  entered  into  the  choice  of  this  function,  it  would  be  applicable  to 
orthotropic  as  well  as  isotropic  materials.  Thus  the  complex  functions  Z^  for  the 
Eq.  (21)  is 


fl<V 

j  “  1.2 

For  a  crack  with  traction  free  surfaces,  the  boundary  conditions  are 

for  x  <  a  and  y  =  0. 

From  Eq.  (13),  note  that  on  y  =  0 


(22) 


T  =0, 

xy 

a  =0, 

y 


(23) 


£i  h  *• 

And  from  Eq.  (22) 

zu  =  zi2  ony  =  °- 

Therefore  the  first  condition  of  Eq.  (23)  is  automatically  satisfied  from  Eq.  (21); 


a 


=  0 


=  (a2  +(32)  Re  (Zn), 


a  =Re(Zn)  =  Re(Z12)  =  Re(Zli). 

y  y  =  0  J 
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In  order  that  the  second  condition  of  Eq.  (23)  be  satisfied, 

=  0, 


y  =  0 
|x|  <  a 


-  Re<V  - Re  n$r 


y  =  0 
lxi<  a 


(24) 


where  H(x)  is  real.  Equation  (24) demands  that 


Im  i1  (S.) 


=  0, 


y  =  0 
(x|  <  a 


j  -  1.2 


(25) 


It  should  also  be  noted  that  f^  (£^)  is  prescribed  by  the  external  loads  far  away 
from  the  crack.  For  the  case  of  a  pure  tensile  traction  o  applied  at  infinity, 


(26) 

It  can  be  seen  that  the  conditions  specified  by  Eq.  (25)  are  satisfied  by  Eq,  (26).  Sub¬ 
stitution  of  Eq.  (26)  and  (22)  into  Eq.  (21)  gives  the  stress  distribution  at'  any  point  in 
the  plate . 

However, if  examination  of  the  stress  intensity  is  limited  to  the  neighborhood  of 
the  tips  of  the  crack,  certain  simplification  is  possible.  Due  to  symmetry,  only  the 
stress  intensity  around  one  crack  tip  need  be  examined. 

For  the  stress  near  the  vicinity  of  the  tip  of  the  crack,  the  location  of  the  geo- 

10 

metrical  points  around  the  tip  of  the  crack  can  be  expressed  as  a  +  rex  where 
0<  |6|<  7r .  Recall  from  Eq.  (13)  that 


S1  =  X+  Sly  ' 

=  X  +  S2y  * 

can  ^  consihered  as  complex  variables  of  two  non -orthogonal  coordinates  with 
non -normalized  base  vectors  s ^  and  s^.  They  can  be  defined  in  the  following  manner:* 


£=  x  +  s,y  =  a  +  r  (cos  0  +  s..  sin  e), 

1  1  (27) 

=  x  +  s2y  =  a  +  r  (cos  0  +  s2  sin 0 ). 

Note  that  and  £2  are  not  the  geometrical  points  at  which  the  state  of  stress  is  under 

consideration.  On  the  other  hand,  in  the  special  case  of  isotropy,  a  =  0  and  (3  =  1  lead  to 

10 

s.^  =  s2  =  i,  or  £  =  a  +  r(cos0  +  i  sin0)  =  a  +  re  ,  i.  e.  £  and  the  geometrical  points  co¬ 
incide  , 
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For  the  points  in  the  vicinity  of  the  tip  of  the  crack, 

r  «  a  (28) 

Eq.  (28)  implies  that 

fx  <9  =  fj_  (a)  (29) 

Substitution  of  Eqs.  (27)  and  (29)  into  Eq.  (22)  leads  to 

Z  =  j  =  1, 2  (30) 

lj  V  2r  (cos  0  +  Sj  sin9  ) 

In  the  sequel ,  the  approximation  notation  is  dropped  for  convenience  but  it  should  be 
kept  in  mind  that  the  expressions  are  applicable  only  in  the  vicinity  of  the  crack  tip. 

Substitution  of  Eq.  (30)  into  Eq.  (21)  gives  the  stress  distribution  around  the  tip 
of  a  crack  in  an  infinite  plane  under  tensile  traction  applied  at  infinity  (Fig.  2a)  as: 


o  .wr  to:„hp 

x  /IF  2 


Ret 


VcosS  +  s^  sinO  VcosG  sinG 


B-Jm  (-—J,  -  _ ! 

a  V cos0  +  s1  sin  G  VcosO+s^sinG 


cr/a"  1 

a  - -  -Tj- 

y  /if  2 


Re  ( 


+ 


VcosW  +  s  sin  G  /cosO  +  's^sinO" 


(31) 


+  Im  ( 


a  VcosG  +  s,  sin  8  Vcose  +s0  sin  e 


T  «  <g2.±P.?> 

xy  /IF  2a 


Re( 


VcosQ  +  Sj^  sin  0  VcosO  +  s2  sine 


b.  Case  2:  Skew -Symmetric  Loading 

For  external  loading  skew-symmetric  to  the  crack,  the  Westergaard- 
Irwin  stress  function  can  be  written  as  * 


U2  "  "  IF 


la  Re  _^21  '  ^22_ 


(32) 


Rewriting  Eq.  (15)  using  the  notation  described  following  Eq.  13. 
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Substitution  of  Eq.  (32)  into  Eq.  (5),  using  Eq.  (18)  and  (19)  for  i  =  2  yields 


[<«2  -  32)  Re  <Zn  -  Z22)  -  2a  0  Im  (Z21  +  Zj,)], 
Re  (Z21  -  Z22), 

Re  (Z21+  Z22)  -  (3  Im  (Z21  -  Z22)J  . 

For  a  single  crack,  the  complex  functions  Z2^  are 

Z, 


CTy  ‘  2a 


T  —  —pz — 
xy  2  a 


f2  (i,) 

'2j= 


j  -  1.2 


The  boundary  conditions  are 


ffy  =0’ 


Txy 


for  lxl^.a  and  y  =  0. 


(33) 


(34) 


(35) 


Considerations  similar  to  those  in  the  symmetric  case  can  be  used.  Equation  (34) 
automatically  satisfies  the  first  condition  of  Eq.  (35),  and  for  the  second  condition  to  hold 


Im  f2  ( £ .) 


y  =  0 
|x!<  a 


=  0 


j  -  1.2 


(36) 


For  a  plate  subjected  to  pure  shearing  traction  at  infinity 


f2<9=T£r  (37) 

Thus  the  approximate  expressions  for  the  complex  functions  defined  by  Eq.  (34)  in  the 
vicinity  of  the  crack  tip  are 


t  /“a 


V"2r*  Vcos  6+s.  sin  0  ’ 
1 


j  =  1,  2  (38) 
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Substitution  of  Eq.  (38)  into  Eq.  (33)  gives  the  stress  distribution  around  the  tip  of  a 
crack  in  an  infinite  plane  under  shearing  traction  at  infinity  (Fig.  2b)  as: 


_  Tyra 

ax  /IF 


2a 


Re( 


Vcos  8  +  s  ^sin  0 


V cos  0  rs^sin  6 


_  t /a~ 
xy“  /IF 


T  = 


+  j3  Im  (■■■t-  .-1' . . ;•  ■  +  —  -y— -  -T .  ) 

Vcos  0+ s  ^lnQ  Vcos  0  +  s„sin  8 


tV2~  1 
CTy  /IF  2a 


Re( 


Vcos  0  +  s^sin  0  Vcos  0  +  s^sin  0 


(39) 


4-  Re(-==== -  +-===4=  _== 

Vcos  0  +  s,  sin  0  Vcos0+  s„sin0 


-  JL  Im  ( 

2 a  V  cos  0  + 


“tT  ) 


V cos  0  +  s^sin  0  Vcos  0  +  s^sin  0 


Since  Eq.  (31)  and  (39)  describe  the  stress  distribution  around  the  tip  of  the  crack  in 
plates  of  identical  geometry  (in  each  case,  a  single  crack  is  present),  the  two  solutions 
can  be  superimposed  to  describe  any  arbitrary  place  traction.  Then  for  the  loading 
condition  as  shown  in  Fig.  2c,  the  crack  tip  stress  distribution  obtained  by  the  Wester- 
gaard -Irwin  method  is 


a  fe*  +£  )  L  Re  (  +  ■  ■  ■  1  )  -  p  lm  (  .  .  .  1  — -  +  .  1  l] 

x  V2r  2a  vcos  0  +  s^sin  6  Vcos  0  +  s2sin  0  Vcos  0  +  s^sine  Vcos  6  +  s2sin  6  J 

-ifEJi  [r „2-p2)  Re  (  •  ■  -  1 - —  -  1 - —\-2nft  Imf- —  1  -  +  1  — — ."3 

V  2r  |_  Vcoa  0  +  s^sln  0  Vcos  6  +  s2sin  0  V  cos  6  +  s^sin  0  Vcos  8  +  SjSin  8  _J 

a  =  _1_  r,  Re  t _ 1 _  +  _ 1 _ 1  +  p  Im  ( _ 1 _  _  _ 1 _ J 

y  2 a  L  Vcos  8  +  SjSin  9  Vcos  6  +  SjSin  9  Vcos  6  +  SjSln  0  Vcos  S  +  s2sin0_| 

tVT  1  Fr  l  .  i  ^1 

■/2F  2 a  [_  Vcos  6  +  Sj^sinB  Vcos  9  +  s2  sin  ej 


(40) 


_  -jT?  (a 
*y  VIET 


2+b2)  rnr,  i _ l___i 

2 a  _  Vcos  6  +  SjSin  8  Vcos  0  +  s2sin  6  _ 


V2F  2a  ^Vcos  0  +  SjSin  9  +  7cosTT+" sin  6  ^  ^  110  ^V cos  9  +’s^sInF  Vcos  0  +  s^sin  0 
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2.  Solution  by  Lekhnitskii's  Method  jVjfel  * 

Returning  to  Eq.  (11) 

u  =  zl(£l)+z2  (£l)  +  Z3  (S2)  +  Z4  (z2)  (11) 

where  t2  are  defined  by  Eq.  (13). 

Corresponding  to  the  stress  function  U  =  Re  [V  ^  (£)  -  ^((£)j  in  the  isotropic  case,  the 
stress  function  for  Eq.  (11)  is 

(41) 


(42) 


(43) 


Lekhnitskii  gives  the  stress  functions  (^)2(£2),  ^3^3)  for  an  anistropic  space 

containing  an  elliptical  cylindrical  cavity  under  arbitrary  external  tractions .  The 
stress  functions  can  be  modified  for  a  plane  containing  an  elliptical  hole.  Performing 
the  necessary  operations,  the  stress  functions  for  the  plane  problem  are: 


U  =  2Re  [Z1(t1)+Z2  (£2)J. 

The  following  additional  notations  are  introduced: 


Substitution  of  (41)  into  (5)  using  (42)  gives 

ox  =  2  Re  (£l)  +  s2  (j)2  (£2)J. 

Oy  =  2  Re  j(p  4(£l)  +  ^2  (£2)^j> 

Txy=  '  2  Re(|  S1  ^1  ^1^  +  S2^2  ^2) _  ' 


00 

3W  =  Vn  +  Sj-hj  2=i  (bm  -  s2  am)  r,{m. 


co 


(44) 


Q2(£2)  a2  ^nr>2  '  s1  -  s2  ^=1  ^bm  '  S1  am  ^  7,2  > 

(^3(t3)  =  0. 


■ - - — - -  V 

*  In  the  following  brief  summary  of  Lekhnitskii's  method,  most  of  the  notations  of  the 
original  author  are  retained.  However, some  of  the  notations  are  changed  in  order  to  re¬ 
late  the  corresponding  terms  which  appeared  previously  in  the  analysis  obtained  by  the 
Westergaard -Irwin  type  of  solution. 
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where 


,  ,VT2 - 2 - rj 

£k  +  £k  -  a  -  sk  b 
_____ 

k 


k  =  1,2 


(45) 


and  a  is  the  major  axis  and  b  is  the  minor  axis  of  the  ellipse. 

For  a  crack  of  length  2a  instead  of  an  elliptical  hole,  the  r)k  in  the  stress  functions 
and  (j)2  ta^e  on  the  special  values  for  b  =  0,  i.e. 


r~2 — t 

£k+  *k  -a 


k  =  1,2 


(46) 


Using  the  notations  defined  by  Eq.  (42)and  letting  b  =  0,  Eq.  (44)  becomes 

^1  ^  =  V  2  1  2  A1  ‘  s  x  s2  *ml  m(Bm  ‘  S2  am>  *1 
*1  "  a 


(47) 


-*-/  l  l  00 

Q2  (t2)  =  7—2 — jr  A2  +  T-ri-  £ml  m(Em  -  s±  am)  n2 

£2  *  a  1  z 


-m 


For  the  loading  condition  as  shown  in  Fig.  2a  Lekhnitskii*  indicated  that  the  constants  in 
Eq.  (47)  are 


a  =  -  acr 
al  2  ’ 

F1-0, 

am=Fm  =  0m^2’ 
Ai  =  a2  =  °. 


For  the  loading  condition  as  shown  in  Fig.  2b,  the  constants  are 


a!  =  °’ 

K  -  aT 

bl  2  ’ 

a  =  b  =0  for  m  >  2  , 
mm  — 


a1  =  a2  =  0. 


This  is  obtained  by  reducing  the  results  of  Ref.  [9]  pp.  129-131  to  a  plane  problem 
and  taking  the  special  case  of  a  degenerate  ellipse  (i.e.  the  minor  axis  b  =  0).  Similar 
results  can  be  obtained  from  Savin's  work  Ref.  [4j  pp.  157-162. 
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For  the  loading  condition  as  shown  in  Fig.  2c,  superimposing  the  two  results  gives  the 

constants  in  Eq.  (47): 

—  a a 

al T  ’ 

T~ 

Hl~  T  ’ 

a  =  B  =0  for  m  >  2, 
mm  —  ' 

Ai  =  A2  =  °. 

Substitution  of  Eqs.  (47)  and  (48)  into  Eq.  (43)  gives  the  stress  distribution  in  the  plate 
under  the  loading  conditions  as  indicated  in  Fig.  2c: 


cr  =  Re 
x 


a  =  Re 

y 


[<V  s2> 


2  2 
-S1  (t  +  s2p  )  +  s2  (t  +  s x cr ) 


r  -  (t  +  s2  o) 


l!V  S2> 


UT/T?' 


’’2 V£22  -  a2  J 
(T  +  S1<T  ) 

-  a2  -J] 


(49) 


t  =  -  Re 
xy 


-s1(T  +  s2a)  s2(t  +  s1ct) 


l<8r  s2)t 


n2V£-2  -  a^ 


To  consider  the  stress  intensity  at  the  vicinity  of  the  tip  of  the  crack  where  r  «  a,  the 
same  simplifications  as  discussed  in  section  II  B  1  can  be  made.  Using  Eq.  (27) 


rji  =  r)2  -  1  when  r  «  a. 

Then  the  crack  tip  stress  distribution  obtained  by  Lekhnitskii's  method  is  * 

S2V  +  sx  o) 

V cos  d  +  s2sin  0 


VT  „  1 

o  = -==.  Re< - - 

x  V2F  (s1-  s2) 


-s1  (t+  s2  a) 
/cos  6  +  SjSin  0 


+ 


X 


a  Re< 

y  rn 


(Si  -  s2) 


‘  -(t  +  s2  a) 

Vcos  0  +  SiSin  0 


+ 


(T  +  B1  CT) 

V  cos  0  +  s2sin8  ^ 


(50) 


T  =—  Re( 
xy  /2r 


(sr  s2> 


sx  (t  +  s2a  ) 
Vcos  9  +  SiSin  0 


s2(r  +  Sla) 


Vcos  9  +  s2  sin  6'  | 


G.  C.  Sih  and  P.C.  Paris  [id)  discussed  a  similar  problem  solved  by  Lekhnitskii.  How¬ 
ever  the  result  given  by  Sih  and  Paris  is  contradictory  to  that  given  here  in  the  respect 
that  the  roles  of  a  and  t  are  reversed  in  addition  to  some  sign  discrepancies. 
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3.  Comparison  of  the  Two  Solutions 

It  should  be  noted  that  both  Eq.  (40)  and  Eq.  (50)  describe  the  stress  dis¬ 
tribution  around  the  tip  of  a  crack  in  an  infinite  orthotropic  plate.  From  physical  con¬ 
sideration,  the  stresses  in  Eq.  (40)  and  Eq.  (50)  are  necessarily  real.  This  indeed  is 
the  case  since  the  real  and  imaginary  parts  of  the  complex  variables  in  the  expressions 
are  real  quantities. 

Further,  both  of  the  expressions  are  simplified  by  neglecting  r  terms  of  order 
-1/2 

higher  than  r  .  Therefore  from  the  uniqueness  theorem  in  the  theory  of  elasticity, 
Eq.  (40)  and  Eq.  (50)  should  be  identical.  Although  little  similarity  can  be  noted  in 
their  present  forms,  they  can  be  compared  by  expressing  the  real  and  imaginary  parts 
differently  as  follows: 

1 _  =  _ 1 _ 

Vcos  9  +  s  sin  6  V(cos  0  +a  sin  0)  +  i  (ft  sin  6) 

1  .  <51> 

. . . : . . :  (  cos  *  *  sin  -5—  ) 

(cos  0  +a  sin  0)2+  ( 0  sin  0)2 


where 


sin<)>1 
cos  (J>1 


(3  sin0 

- ^ > 

(cos0  +a  sin0)  +  <J3  sin0) 

cos  0+a  sin  0 

V  2  _  . . ~2  ’ 

(cos  0  +a  sin  0)  +  ((3  sin  0) 


(52) 


similarly 


where 


<t>^ 


cos  0  +  s„  sin  0 


- : -  (cos  ^-i  sin  T  ) 

cos  0  -  a  sin  0  )  +  (0  sin  0)  (53) 


sin<j>2 

cos<j>2 


_ [ 3  sin  0 _ 

V - 2 - T  ’ 

(cos  0  -a  sin  0  )  +  (j3  sin  0) 

cos  0  -  a  sin  0 _ 

- 2 - j  • 

(cos  0  -  a  sin  0  )  +  ((3  sin  0) 


(54) 
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Using  Eqs.  (51)  to  (54),  it  can  be  shown  that  Eq.  (40)  and  Eq.  (50)  are  identical  and  they 
can  both  be  expressed  as 


x  /F  2a 


.  Tv^~  /  1 

Vlr  2a 


♦  l  4>2 

a  cos  -j-  +  |B  sin  -y 


v: 


1 - 2 

]_'(cos0+a  sinG)  +  ((3sin0) 

4>i  i 


+ 


^2  ^2 
acos-y-  -  j3sin-y 


*  * 

(cosG  -a  sinG)  +(j3  sinG)  __ 

'9  2  Yi  ti  99  4*0-1 

(a- (3  )cos  -j—  +  2a  (3sin-y-  (a  -(3  )cos  -y  +  2a (3sin  -y- 


V  2  2" 

(cosG+a  sinG)  +  ((3 sinG) 


/  'Ar  rt 

v(cosG -a  sinG)^ +(3sin0)z  J 


cr/T,  1  . 


xVa".  1  . 


4>i  4=! 

a  cos  -75 —  (3  sin 


V7 


2  2 

(cosG+a  sinG)  +(J3  sinG) 


4>i 


cos 


4*2  4*2 

a  cos  ~2 — +j3  sin  -y- 

(cosG-a  sinG)  +(j3sinG)  _| 

4>i  n  (55) 

COS 


v - 2 - 

(cosG+a!  sinG)  +  (/3  sinG) 


V  2  2 

(cosG  -a  sinG)  +  (/3sinG) 


xy  V  2r  2  a 


cos 


"2“ 


cos 


r  “  A  p  /  ft  . . 

v(cosG+a  sinG;  +(j3  sinG)  v(cosG  -a  sinG)  +(|3  sinG; 


V  2r  2a 


acos 


1  , 


(3  sin  -5— 


4*2  ^2 
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where  4^  andcj>2  are  defined  by  Eq.  (52)  and  (54). 


In  general,  arbitrary  external  plane  tractions  can  be  resolved  in  terms  of  c r^, 
a2  and  t  as  shown  in  Fig.  lb.  However,  it  can  be  shown  [4]  that  the  stress  distribution 
in  an  anisotropic  plate  is  not  influenced  by  a  straight  crack  in  the  direction  of  the  tensile 
stress.  In  other  words,  the  component  of  tensile  traction  in  the  direction  of  a  crack 
does  not  contribute  to  the  stress  singularity  at  the  crack  tip  and  its  effect  on  the  stability 
of  the  crack  is  small  and  can  be  neglected.  Consequently  Fig.  2  c,  in  fact,  represents 
the  general  case  of  plane  loading  in  the  absence  of  bending  and  twisting.  For  a  thin  ortho - 
tropic  plate  containing  a  single  crack  in  the  direction  of  one  of  its  two  axes  of  elastic 
symmetry  under  an  arbitrary  external  tractions,  Eqs.  (55)  gives  the  stress  distribution 
around  the  tip  of  the  crack. 
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C.  Comparison  of  the  Crack  Tip  Stress  Distribution  in  an  Orthotropic  Plate  to 
that  in  an  Isotropic  Plate. 

For  an  isotropic  plate  containing  a  single  crack  under  the  loading  condition 

1/2 

as  shown  in  Fig.  2c,  and  after  neglecting  r  terms  of  order  high  than  r  ,  the  crack  tip 
stress  distribution  can  be  expressed  as 
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where  the  stress  intensity  factors  are 


k^  =  oVT  (symmetric)  , 

k2  =  t7a~  (skew-symmetric). 

Comparison  of  Eq.  (55)  to  Eq.  (56),  leads  to  the  following  observations.  The 

most  obvious  difference  between  the  state  of  stress  in  the  vicinity  of  the  tip  of  a  crack  in 

an  orthotropic  plate  and  that  in  an  isotropic  plate  is  that  the  former  is  dependent  on  the 

elastic  constants  of  the  material*  while  the  latter  is  not.  However,  there  are  also 

several  similarities.  First,  the  stress  singularity  at  the  tips  of  the  crack  (r  approach¬ 
'll 

ing  zero)  for  both  cases  is  of  the  order  r  '  .  Second,  a  parameter  characterizing 
the  effect  of  the  size  of  the  crack  and  the  magnitude  of  the  external  stress  for  the  ortho - 


The  quantities  a,  ^ are  functions  of  Ex>  Ey,  vx>  and  Gxy. 
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tropic  plate  can  be  defined  in  a  manner  similar  to  the  isotropic  case,  namely, 


=  crZT, 
k2  =  tVT. 


(57) 


The  stress  intensity  factors,  k^  and  k2  therefore  can  be  used  as  parameters  to  measure 
the  strength  of  a  cracked  orthotropic  plate.  The  crack  extension  force,  and  2 
obtained  from  energy  considerations  by  Irwin  [2]  are  generally  considered  to  be  more 
fundamental  parameters  related  to  crack  extension.  However, in  order  to  obtain  ^ 
and  tj  2  i  a  knowledge  of  the  direction  of  the  crack  growth  is  necessary.  In  an  iso¬ 
tropic  material  a  generally  accepted  hypothesis  which  postulates  that  crack  growth  is 
perpendicular  to  the  direction  of  greatest  tension  has  been  verified  experimentally. 
However,  no  such  hypothesis  is  possible  for  orthotropic  materials;  since  the  direction  of 
crack  growth  depends  on  the  strength  of  the  material  as  well  as  the  state  of  stress.  Due 
to  this  difficulty,  the  crack  instability  of  orthotropic  plates  will  be  examined  using  the 
stress  intensity  factors  as  criteria. 


Defining  k^c  and  k2c  as  the  critical  magnitude  of  the  stress  intensity  factors  for 
unstable  crack  extension  under  simple  tension  and  under  pure  shear  respectively,  the 
combined  influence  of  stress -intensity  factors  can  be  expressed  in  the  following  form: 


k.  m 

<kT>  +  < 
lc 


1. 


(58) 


Obviously,  in  order  to  express  k^and  k2cin  a  functional  relationship,  a  fracture  criterion 
relating  the  influence  of  symmetric  and  skew -symmetric  loading  for  orthotropic  material 
is  required.  At  present,  not  only  is  a  fracture  criterion  of  this  type  unavailable  for 
orthotropic  materials,  but  a  workable  criterion  has  yet  to  be  formulated  and  validified  for 
homogeneous  isotropic  materials.  Consequently  the  exponents  m  and  n  in  Eq.(58)  must 
be  determined  experimentally.  From  the  experimental  results,  a  suitable  fracture 
criterion  for  orthotropic  material  may  become  evident. 
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III.  EXPERIMENTAL  PROGRAM 

A.  Introduction  to  the  Experimental  Investigation 

1.  Purpose  of  the  Experimental  Program 

An  experimental  program  was  designed  to  determine  the  validity  of  the 
hypothesis  formulated  from  the  analytical  results.  There  are  two  main  objectives  to 
the  experimental  program  for  orthotropic  plates:  (a)  To  investigate  the  feasibility  of 
using  stress  intensity  factors  and  as  criteria  for  unstable  crack  extension,  and  (b) 
To  investigate  if  a  functional  relation  between  k^  and  k^  exists  and,  if  it  does  exist, 
whether  it  can  be  expressed  in  the  form  of  Eq.  (58). 

The  material  chosen  for  the  experimental  investigation  was  balsa  wood.  This 
choice  was  based  on  two  factors;  first  because  it  was  essentially  orthotropic, and  second 
because  it  was  readily  available  and  suitable  specimens  could  be  easily  prepared  and 
tested.  However,  due  to  the  non -uniformity  of  the  wood*  only  a  qualitative  demonstra¬ 
tion  of  the  fracture  mechanics  to  orthotropic  plates  should  be  expected. 

2.  Experimental  Set-up 

In  these  experiments,  the  desired  loading  can  be  classified  into  three  types: 
(a)  Pure  tension  perpendicular  to  the  crack,  (b)  Pure  shear  and  (c)  Combined  tension  and 
shear.  Two  methods  of  testing  were  used  to  achieve  these  loading  conditions. 

In  the  first  method  of  testing,  the  tension  test  as  shown  in  Fig.  3a,  loading  from 
pure  tension  perpendicular  to  the  crack  to  combined  tension  and  shear  can  be  obtained  by 
varying  the  relative  angle  y  between  the  applied  tension  and  the  line  of  the  crack  **. 

The  specimen  was  gripped  at  opposite  ends,  one  end  being  pinned  to  a  table  and  allowed 
to  rotate.  The  tension  load  was  transferred  to  the  other  grip  through  a  long  cable  to 
minimize  any  bending  moment  induced  by  the  deformation  of  the  specimen.  However, 
pure  shear  could  not  be  obtained  by  this  method.  Consequently  a  second  loading  method 
was  used. 

A  dead-weight  testing  machine  shown  in  Fig.  5  was  constructed  to  load  the  speci¬ 
men  in  pure  shear  as  well  as  in  combined  tension  and  shear.  Its  essential  components 
are  shown  in  the  schematic  drawing  in  Fig.  6.  A  square  balsa  wood  specimen  was 
gripped  on  all  four  sides.  Grip  No.  1  was  pinned  to  the  machine  base  so  that  rotation 
in  the  plane  of  the  plate  was  possible.  Cables  connected  to  the  remaining  three  grips 

*  This  will  be  discussed  in  a  later  section, 

**  Since  it  has  been  noted  in  section  II  B  3  that  tension  in  the  direction  of  the  crack  has 
no  effect  on  the  stress  singularity  at  the  crack  tip,  a  in  Fig .3a  need  not  be  considered 
and  the  loading  condition  can  be  considered  to  be  identical  to  that  represented  in  Fig.  2c. 
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were  arranged  as  shown  in  Fig,  6;  the  three  cables  were  then  connected  together  in  a 
manner  such  that  the  load  on  a  single  center  pan  was  distributed  equally  to  each  of  the 
three  cables.  These  loads  combined  with  the  reaction  force  from  the  pin  applied  a  pure 
shear  to  the  specimen.  The  load  transferring  cables  were  of  sufficient  length  so  that 
the  pure  shear  condition  remained  unaffected  when  the  specimen  deformed  under  load. 
For  tension,  another  cable  was  attached  to  grip  No.  2  (in  Fig.  6)  and  loaded  by  a  sepa¬ 
rate  weight.  The  detailed  manner  in  which  the  cables  were  connected  to  the  grips  is 
shown  in  Fig.  7.  Note  that  the  pin  (A)  was  located  in  the  center  of  the  gripping  area  so 
that  no  moment  was  produced.  The  construction  of  the  grips  is  also  shown  in  Fig.  7. 
Each  grip  was  assembled  from  two  aluminum  pieces  .  The  gripping  surfaces  were  made 
of  pads  of  rubber  which  has  an  elastic  modulus  close  to  that  of  the  balsa  wood  thus  mini¬ 
mizing  the  constraint  of  the  grips .  By  means  of  a  torque  wrench,  a  minimum  tightening 
force  required  to  prevent  the  grips  from  slipping  under  load  was  determined  and  main¬ 
tained  uniformly  from  specimen  to  specimen. 

To  prevent  buckling  in  shear,  the  grips  were  constrained  in  a  plane  by  two  plates 
of  glass.  Plexiglass  plates  were  affixed  to  the  glass  plates  to  further  constrain  the 
specimen  from  buckling.  A  cross-sectional  view  indicating  their  relative  position  is 
shown  in  Fig.  8.  The  clearance  between  the  plexiglass  plates  and  the  specimens  was 
approximately  0.  005  in.  The  friction  between  the  aluminum  grips  and  the  glass  plates 
was  very  small  and  the  displacement  of  the  grips  was  not  significantly  restricted. 

Observation  of  the  crack  was  facilitated  by  placing  a  light  below  the  specimen. 
The  crack  length  was  measured  using  a  traveling  microscope  of  20X  power  to  which  a 
vernier  scale  was  attached  thus  providing  length  measurements  accurate  to  0.01  inch, 
Fig.  5  and  8. 

B.  Testing  Procedure 

The  stress  analysis  presented  in  Part  I  is  applicable  to  a  single  crack  parallel 
to  one  of  the  principle  planes  of  elastic  symmetry  of  a  large  plate.  In  these  experi¬ 
ments,  only  the  case  of  a  crack  parallel  to  the  grain  of  the  balsa  wood  was  investigated 
since  it  is  the  most  frequently  observed  mode  of  cracking  in  orthotropic  plates. 

Four  different  combinations  of  load  were  employed.  Referring  to  Fig.  9,  path 
No.  1  represents  tension  perpendicular  to  the  crack.  Path  No.  2  represents  a  crack 
orientated  at  an  angle  y  to  the  line  of  tension.  Path  No.  3  represents  a  specimen  loaded 
first  to  a  pure  tension  stress,  <j  ,  and  then  loaded  to  fracture  by  increasing  the  shear 
stress.  Path  No.  4  represents  pure  shear  loading.  As  stated  in  the  previous  section, 
two  types  of  testing  methods  were  used:  first  a  tension  test  with  the  crack  orientated  at 
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various  angles  with  respect  to  the  line  of  tension  and  second  a  pure  shear  test  or  com¬ 
bined  tension  and  shear  test  using  four  grips  to  apply  the  desired  loads.  Since  different 
procedures  for  the  two  types  of  testing  are  used,  they  will  be  discussed  separately. 

1.  Tension  Tests  (Path  No.  1  and  2) 

Five  groups  of  specimens  were  tested  in  tension.  The  dimensions  and 
the  values  of  the  angle  y  for  each  group  are  listed  in  the  table  below. 


Group 

y-deg 

Thickness 

Width 

Length 

in 

in 

in 

1 

90° 

1/16 

6 

6 

2 

65° 

1/16 

4-1/2 

6 

3 

50° 

1/16 

4-1/2 

6 

4 

35° 

1/16 

4-1/2 

8-1/2 

5 

20° 

1/16 

4 

14 

The  variations  in  dimensions  of  the  specimens  were  dictated  by  two  considerations:  (a) 
The  longer  lengths  were  provided  for  small  values  of  y  so  as  to  prevent  crack  propa¬ 
gation  into  the  grips,  (b)  The  narrower  widths  of  the  specimens  were  necessitated  by 
the  size  of  the  balsa  wood  plate  available. 

An  initial  crack  was  put  in  the  center  of  the  specimen  which  was  clamped  in  the 
grips  and  loaded  in  tension  by  dead-weights.  The  crack  length  was  measured  using  the 
microscope,  and  the  load  increments  at  which  crack  length  measurements  were  made 
were  large  at  low  loads  and  reduced  accordingly  when  the  crack  length  was  observed  to 
increase.  This  procedure  was  continued  until  fracture  occurred.  After  the  specimen 
was  fractured,  it  was  observed  that  the  crack  usually  did  not  propagate  at  the  original 
angle  y,  but  at  a  slightly  different  angle.  The  actual  angle  of  fracture  cf  the  region 
near  the  center  of  the  specimen  was  measured.  From  this  angle,  the  component  of 
tension  stress  ct  perpendicular  to  the  crack  and  the  shear  stress  component  t  were 
computed.  Typical  results  of  the  experiments  are  shown  in  Fig.  10  and  11  where  the 
observed  crack  lengths  2a  are  plotted  as  abscissa  corresponding  to  the  value  of  a  or  t, 
shown  as  ordinates.  These  results  will  be  discussed  later. 

2.  Combined  Tension  and  Shear  and  Pure  Shear  Tests  (Path  No.  3  and  4) 

In  the  second  series  of  tests,  all  of  the  specimens  had  dimensions  of  6"  x 
6"  x  1/16".  From  the  construction  of  the  testing  machine,  it  was  clear  that  the  side 
grips,  No.  3  and  4  in  Fig.  6,  must  be  in  place  in  order  to  apply  shear  load.  However, 
if  tension  load  was  applied  after  the  side  grips  were  in  place,  these  side  grips 
would  constrain  the  elongation  of  the  specimen.  Consequently  in  the  combined  tension 
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and  shear  tests  which  follow  path  No.  3  a  different  procedure  was  used.  For  each 
group  of  specimens,  the  tension  stress  was  increased  to  the  desired  level  with  the  side- 
grips  absent.  Then  the  side  grips,  No.  3  and  4,  were  put  in  place  and  shear  load  was 
increased  while  the  crack  length  was  measured  at  increments  up  to  fracture.  Four 
such  groups  of  specimens  were  tested  for  constant  tension  stress  of  25,  48.6,  62  and 
71  psi  respectively.  The  typical  relations  of  x  vs  2a  observed  at  constant  cr  are  shown 
in  Fig.  12. 

In  the  pure  shear  test,  a  crack  was  located  in  the  center  of  the  specimen  parallel 
to  the  wood  grain  and  its  length  was  measured  at  appropriate  load  increments.  A 
typical  diagram  of  x  vs  2a  is  presented  in  Fig.  13. 

C.  Interpretation  of  Data  and  Results 

In  the  experimental  investigation,  it  was  observed  that  in  almost  all  cases,  some 
stable  crack  extension  took  place  prior  to  the  sudden  fracture  of  a  specimen.  The 
amount  of  stable  crack  extension,  the  load  level  at  which  the  crack  began  to  extend,  and 
the  rate  of  extension  all  vary  from  specimen  to  specimen.  Slight  variations  in  the 
properties  of  the  wood  are  believed  to  be  the  cause  of  these  differences.  Cracks  in  the 
majority  of  the  specimens  tend  to  start  extending  at  approximately  85%  of  the  load  level 
at  which  rapid  fracture  occurs.  Typical  data  for  specimens  No.  52,  110  and  44  are 
shown  in  Fig.  10,  11  and  13  respectively.  Other  frequently  observed  crack  extension 
patterns  are  described  below. 

A  crack  which  was  located  in  a  low  density  region  of  the  wood  normally  extended 
much  more  than  one  located  in  a  higher  density  region.  Examples  of  each  type  of  be¬ 
havior  are  specimens  No.  105  and  41  shown  in  Fig.  11  and  13.  Sometimes,  a  defect  in 
the  wood  was  located  near  the  tip  of  the  crack.  In  this  case,  the  crack  began  to  extend 
at  a  low  load  level,  but  stopped  after  it  had  propagated,  through  the  defective  region. 

This  crack  remained  at  a  constant  length  until  the  fracture  load  was  approached  when  it 
again  resumed  the  normal  extension  pattern  to  fracture.  Specimen  No.  14  in  Fig.  12  is 
a  typical  example  of  this  pattern  of  crack  extension.  Another  irregularity  observed 
was  an  obstacle  which  appeared  in  the  microscope  as  a  dark  spot  located  in  the  path  of 
the  crack.  An  extraordinarily  high  stress  level  was  required  for  the  crack  to  over¬ 
come  the  obstacle.  On  the  other  hand,  a  very  low  fracture  stress  resulted  when  the 
crack  extended  into  a  large  weak  region  such  as  a  long  weak  grain.  However,  the 
worst  irregularities  were  avoided  by  judicious  placement  of  the  crack  in  a  grossly 
homogeneous  region  by  examining  the  specimen  in  front  of  a  strong  light  source. 

Since  the  purpose  of  the  experiments  was  to  investigate  the  unstable  crack  propa- 
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gation  phenomena,  the  desired  information  was  the  stress  level  and  crack  length  at 
which  the  crack  propagated  rapidly  to  fracture.  Because  various  crack  extension 
patterns  were  observed,  a  suitable  convention  was  required  so  that  the  critical  stress 
level  and  the  critical  crack  length  could  be  defined.  The  stress  at  fracture  was  defined 
as  the  critical  stress  level.  On  diagrams  of  cr  vs  2a  or  t  vs  2a, a  straight  line  was 
drawn  tangent  to  the  last  measured  stable  crack  length  and  its  intercept  with  the  critical 
stress  level  was  defined  as  the  critical  crack  length.  In  effect,  this  definition  of 
critical  crack  length  assumed  that  the  stable  crack  extension  rate  was  unchanged  from 
the  last  measurement  up  to  the  point  of  fracture.  This  assumption  was  reasonable  be¬ 
cause  the  last  incremental  increase  of  load  was  usually  very  small  as  can  be  observed 
from  the  graphs. 

As  stated  earlier  the  fracture  angles  in  the  tension  tests  often  differed  slightly 
from  the  original  crack  orientation  angle.  This  variation  was  observed  in  each  group 
of  tests.  To  account  for  this,  the  arithmetic  mean  of  the  final  fracture  angles  of  all 
specimens  tested  in  each  group  was  computed  and  denoted  by  Tmean-  The  results  of 
each  group  of  specimens  are  shown  on  diagrams  of  log  o  vs  log  a  or  log  t  vs  log  a. 
Referring  to  the  loading  paths  in  Fig.  9,  the  results  of  tests  which  follow  path  No.  1  is 
presented  in  Fig.  14,  and  the  results  of  the  tests  which  follow  path  No.  4  is  presented 
in  Fig.  15.  The  results  of  path  No.  2  at  y  equal  to  64°,  50°,  35.  5°  and  22°  are 
presented  in  Fig.  16a  and  b,  17a  and  b,  18a  and  b,  19a  and  b,  respectively.  Finally, 
the  results  of  path  No.  3  at  constant  tension  stress  ctq  equal  to  25,  48.6,  62,  and  71 
psi  are  presented  in  Fig.  20,  21,  22, and  23  respectively. 

D.  Discussion  of  Results 

In  the  analysis  of  stress  in  the  vicinity  of  the  tips  of  a  crack,  it  was  suggested 
that  the  stress  intensity  factors  k1  and  k2  could  be  used  as  parameters  to  measure  the 
strength  of  an  orthotropic  plate  containing  a  crack.  If  this  is  feasible  (as  it  is  in  the 
isotropic  case)  the  critical  stress  intensity  factors  k^c  and  k2c  could  be  constants  for 
one  orthotropic  material.  Equation  (57)  can  be  written  in  the  following  form: 

log  k  =  logo  +  -i-  log  a, 

1C  (59) 

logk2c  =  log  t  H  log  a. 

If  k  and  k2c  are  constants,  then  Eq.(59)  suggests  a  linear  relation  between  logo  and 
log  a  or  log  t  and  log  a  with  slope  equal  to  -0.5  on  log-log  plots. 

In  order  to  examine  the  hypothesis  that  the  experimental  results  represent  a 
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straight  line  of  slope  equal  to  -0.  50,  the  data  obtained  from  the  tests  which  follow  path 
No.  1  are  presented  in  Fig.  14.  The  a  and  a  of  all  of  the  experimental  data 
were  computed  and  are  indicated  in  Fig.  14  by  a  cross  (x).  Next  the  mean  of  the  data 
in  the  region  a  <  a  and  similarly  the  mean  of  the  data  in  the  region  a  >  a_ 
(indicated  by  two  solid  dots)  were  located.  If  a  straight  line  passes  through  the  three 
points,  it  is  reasonable  to  conclude  that  these  data  may  be  represented  by  a  straight 
line.  In  Fig.  14  it  can  be  seen  that  the  three  points  fall  on  a  straight  line  the  slope  of 
which  was  determined  graphically  to  be  -0.47.  However,  since  a  slope  of  -0.  5  was 
predicted  by  the  analysis,  a  line  of  this  slope  was  drawn  through  the  mean.  A  band 
of  +  10  per  cent  about  this  line  included  almost  all  of  the  experimental  data.  Within 
the  realm  of  experimental  scatter,  the  critical  stress  intensity  factor  k^  =  a  -fH  = 
constant  provided  a  good  description  of  the  data  for  orthotropic  plates  tested  under 
loading  symmetrical  to  a  crack. 

In  the  same  manner,  the  results  of  the  pure  shear  experiments  which  are  pre¬ 
sented  in  Fig.  15  were  analyzed,  and  these  data  also  were  well  represented  by  a 
straight  line  with  a  slope  equal  to  -0.48.  Again  a  band  of  +  10  per  cent  about  the  line 
of  -0.  5  slope  through  the  mean  included  most  of  the  data.  Consequently,  it  was  con¬ 
cluded  that  the  critical  stress  intensity  factor  k2c  =  xV~a  =  constant  is  likewise  applic¬ 
able  to  orthotropic  plates  tested  under  loading  skew  symmetrical  to  a  crack. 

In  the  tests  which  followed  path  No.  2  where  a  tension  load  was  applied  at  an 
angle  y  to  the  crack,  the  tension  component  a  perpendicular  to  the  crack  and  the  shear 
component  t  were  calculated  from  the  actual  fracture  angle  of  each  specimen.  If  the 
fracture  angles  of  the  specimens  in  each  group  were  the  same  then  the  ratio  of  ct/t 
would  have  been  constant.  But,  since  these  fracture  angles  varied  slightly  from  speci¬ 
men  to  specimen  in  the  same  group,  it  is  evident  from  Mohr's  circle  that  the  ratio  of 
o/t  is  not  constant  and  the  data  would  be  distributed  differently  in  the  graphs  of  log  a  vs 
log  a  than  the  data  represented  in  the  graphs  log  r  vs  log  a.  Thus  it  is  necessary  that 
the  data  be  represented  in  both  of  these  graphs, and  they  are  presented  in  Figs.  16  to  19. 
The  measured  mean  fracture  angles,  the  number  of  specimens  and  the  slopes  of  straight 
lines  estimated  from  these  data  are  summarized  in  the  following  table. 
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Original  Angle 

Mean  Angle 
of  F  racture 

.  No.  of  Specimens 
Tested 

k2=  (a)  a 

k2=  (a)  t 

y 

.  n 

y 

mean 

N 

h 

2 

65 

64° 

8 

-.65 

-.57 

50° 

50° 

7 

-.40 

-.43 

35° 

35.5° 

12 

-.55 

-.53 

20° 

22° 

10 

-.40 

-.46 

Because  of  the  increased  scatter  and  the  relatively  few  data  points,  the  straight  lines 
that  were  estimated  from  the  data  exhibited  slopes  h  and  2  which  differed  measurably 
from  -0.  50.  However  further  examination  of  the  data  revealed  that  the  slopes  of  those 
lines  which  were  determined  from  the  larger  number  of  data  points  were  closer  to  -0.50. 
Furthermore,  the  average  values  of  h  and  2  of  the  four  groups  of  tests  were  both  very 
nearly  equal  to  -0.50.  Therefore  it  appeared  that  as  the  number  of  data  points  were 
increased,  the  slope  of  the  straight  line  determined  by  these  data  would  closely  approach 
-0.  50.  Within  the  realm  of  the  accuracy  of  the  experiments,  straight  lines  of  slope 
equal  to  -0.  50  were  drawn  through  the  point  determined  by  crmean  and  amgan  (indicated 
by  a  cross)  in  each  figure.  A  band  enclosing  a  region  +  10  per  cent  about  the  mean 
line  included  a  majority  of  the  data.  Thus  these  data  provide  additional  evidence  that 
a  straight  line  of  slope  -0.  50  was  a  good  description  of  the  relation  between  log  o  and 
log  a  and  between  log  t  and  log  a  as  suggested  by  Eq.  (59).  From  these  data,  two  im¬ 
portant  conclusions  emerge  regarding  an  orthotropic  plate  which  contains  a  crack  and  is 
tested  under  combined  tension  a  and  shear  t.  First,  the  stress  intensity  factors 

=  crVT  and  k^  =  TV~a  provide  a  good  description  of  the  relation  between  a  and  a,  and 
t  and  a  for  rapid  crack  extension.  Second,  various  combinations  of  ^  and  k2  may  be 
considered  as  critical  combinations  which  produce  unstable  crack  extension  to  fracture. 
Thus  these  data  suggest  that  a  function  f(k^,  k^)  =  ^cr|tica^  exists,  which  describes  the 
conditions  for  unstable  fracture . 

The  results  of  the  combined  tension  and  shear  experiments  which  followed  path 
No.  3  are  presented  in  Fig.  20  to  23.  For  combined  tension  and  shear,  the  critical 
stress  intensity  factor  k^c  =  oV a  for  rapid  crack  extension  had  been  shown  to  be  a  con¬ 
stant.  Thus  for  experiments  described  by  path  No.  3,  if  o  is  constant,  then  k^  cannot 
be  constant.  Therefore  the  data  of  the  tests  which  followed  path  No.  3  representing 
log  t  as  a  function  of  log  a  do  not  produce  a  straight  line,  but  rather,  a  curve  repre¬ 
senting  f  (kr  k2)=  Critical. 

To  examine  the  function  f  (k^  k^)  =  fcritjcaj>  that  is  the  interaction  between  ^ 
and  k2  under  combined  loadings,  the  values  of  k^  and  k2  are  needed.  Since  the  relations 
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=  cr-/"a  and/or  =  tVT  were  applicable  for  all  loading  conditions,  tension,  pure 
shear  and  combinations  of  tension  and  shear,  the  values  of  k.^  and  k^  for  all  of  the  ex¬ 
perimental  data  were  computed.  These  values  are  shown  in  Fig.  24  where  k^  is  the 
ordinate  and  k^  is  the  abscissa,  The  data  in  Fig.  24  define  the  functional  relationship 
between  k^  and  k^  corresponding  to  unstable  crack  extension.  A  curve  is  drawn 
through  the  points  to  represent  the  function  f(k^,  k^)  =  fc>  with  end  points  for  tension 
f(k1>  0)  =  klc  and  for  pure  shear  f(0,  k^)  =k2c>  A  band  which  represents  the  maximum 
deviation  from  the  mean  curve  resulting  simultaneously  from  k^  +  10  per  cent  and 
k^  +  10  per  cent  is  also  indicated.  For  the  purpose  of  illustration  this  average  curve 
and  the  deviation  band  are  replotted  in  Fig.  20  to  23. 

To  obtain  the  function  f(k..,  k2)  =  f,,  a  dimensionless  representation  was  used. 
Using  the  curve  in  Fig.  24  which  represents  the  average  of  the  experimental  data,  a 
graph  of  k1Alc  vs  k2/k2c  is  plotted  in  Fig.  25.  The  curve  relating  k^  and  k2  was 
assumed  in  the  form  of  Eq.  (58): 


The  exponents  m  and  n  were  determined  by  choosing  two  points  on  the  curve  in  Fig.  25. 
Substituting  these  into  Eq.  (62)  led  to  two  simultaneous  equations  from  which  m  and  n 
were  obtained.  This  procedure  was  repeated  several  times  using  different  pairs  of 
points  on  this  curve  until  an  average  value  of  m  and  n  was  obtained.  The  values  of  m 
and  n  which  were  found  to  fit  the  experimental  data  were 

m  =  1. 03 , 
n  =  1.88. 


The  important  features  of  this  function  with  the  above  values  of  m  and  n  are  that  the 
exponents  m  and  n  are  not  equal, and  the  function  is  defined  only  in  the  domain  k^  >  0 
and  k2  >  0  since  m  and  n  are  not  integers.  Referring  to  Fig.  26,  it  can  be  seen  that 
f(k1(  k2)  =  f£  is  necessarily  even  in  k2>  i.  e.  f^k^  ~  k2)  =  f£  (k^,  k2)  because  the 
direction  of  shear  is  arbitrary.  Furthermore  under  combined  compression  and  shear, 
the  curve  f(k^,  k2)  =  f£  is  likely  to  extend  along  the  broken  line  shown  in  Fig.  26.  To 
meet  these  physical  requirements,  the  exponents  m  =  1,  n  =  2  were  chosen.  This 
gave  a  function  that  is  even  in  k2  and  defined  for  k^  >  0  as  well  as  k^  <  0.  Then  Eq.  (58) 


becomes 


(60) 
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Equation  (60),  which  is  shown  by  the  broken  curves  in  Fig.  24  and  25  defines  a  curve 
that  is  not  significantly  different  from  the  average  experimental  curve.  This  equation 
fits  the  data  excellently  and  satisfies  the  above  mentioned  physical  considerations. 
Consequently,  Eq.  (60)  was  used  to  represent  the  function  f(k^,  k£)  =  f  for  the  balsa 
wood  sheets  tested. 

E .  Discussion  of  Experimental  Accuracy 

The  major  sources  of  error  in  the  experimental  investigation  were:  1.  The  non- 
uniformity  of  the  balsa  wood  used  for  specimens.  2.  The  accuracy  of  the  crack  mea¬ 
surements.  3.  The  constraint  introduced,  by  the  side  grips  in  the  combined  tension  and 
shear  tests.  4.  Buckling  of  the  specimen  under  shear.  5.  Absolute  similitude  was  not 
maintained  between  the  theoretical  model  in  the  analysis  and  the  model  used  for  experi¬ 
mental  investigation.  These  sources  of  error  will  be  considered  separately. 

1.  The  non-uniformity  of  the  balsa  wood  used  for  specimens.  The  non - 
uniformity  of  wood  can  be  classified  into  two  types,  first,  localized  non -homogeneity 
of  wood  which  probably  caused  the  different  crack  extension  patterns  described  in 
section  III  C  and  second,  discrete  macroscopic  regions  of  different  densities  within  a 
specimen  as  shown  by  the  light  and  dark  bands  in  the  picture  of  a  typical  specimen  in 
Fig.  4a. 

It  was  found  that  by  examining  the  specimen  in  front  of  a  strong  light  source, 
the  non -homogeneous  regions  could  be  observed  and  the  density  of  the  wood  could  be 
compared.  Using  this  technique,  balsa  wood  sheets  which  showed  approximately  the 
same  uniform  density  in  the  prospective  crack  regions  were  chosen  for  testing.  This 
method  of  choosing  specimens  eliminated  the  major  variations  in  wood.  An  investi¬ 
gation  was  designed  to  determine  the  influence  of  the  specific  gravity  on  the  fracture 
strength  of  the  balsa  wood  specimens . 

Since  discrete  regions  of  different  densities  exist  within  a  specimen,  it  is  the 
specific  gravity  in  the  immediate  vicinity  of  the  crack  rather  than  the  average  specific 
gravity  of  the  entire  specimen  that  should  be  examined.  After  a  specimen  had  been 
tested,  strips  of  wood  approximately  1/4"  wide  on  each  side  of  the  fracture  surfaces 
were  cut  out  and  the  specific  gravity  of  these  two  strips  were  measured.  This  method 
of  specific  gravity  measurement  was  applied  to  specimens  which  were  tested  in  tension 
perpendicular  to  the  crack  and  those  tested  in  pure  shear.  The  deviation  from  the 
mean  specific  gravity  was  denoted  by  AS.G.  The  deviation  of  the  critical  stress  in¬ 
tensity  factor  of  each  specimen  from  the  respective  means  are  denoted  by  Akc-  A 
plot  of  Akcvs  AS.G.  is  shown  in  Fig.  27  where  the  solid  points  and  open  points  represent 
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data  obtained  from  the  tension  test  and  pure  shear  test  respectively.  In  this  figure  the 
trend  indicates  that  the  variation  in  the  critical  stress  intensity  factors  was  directly  re¬ 
lated  to  the  specific  gravity.  This  trend  is  well  defined  indicating  that  the  variation 
in  specific  gravities  was  a  major  factor  which'  contributed  to  the  scatter  of  critical 
stress  intensity  factors.  However,  no  attempt  was  made  to  correct  the  experimental 
data  because  the  scatter  in  Fig.  27  makes  clear  that  any  feasible  method  of  measuring 
specific  gravity  would  introduce  an  average  value  over  a  finite  volume  of  wood.  This 
would  improve  the  appearance  of  the  data  but  it  also  would  introduce  another  variable, 
namely  the  volume  of  the  strip  used  in  specific  gravity  measurements.  Since  the 
specific  gravity  at  the  crack  tip  could  not  be  measured,  the  data  was  used  without  an 
arbitrary  correction.  In  its  present  form,  the  data  showed  the  amount  of  scatter  that 
must  be  expected  in  balsa  wood.  Had  a  more  uniform  orthotropic  material  been  chosen 
for  testing,  scatter  of  the  type  shown  in  Fig.  24  undoubtedly  would  be  reduced. 

2.  The  accuracy  of  the  crack  measurements.  The  measurements  are 
made  through  a  single  microscope.  The  location  of  one  tip  of  the  crack  was  first  deter¬ 
mined  and  the  microscope  was  then  moved  to  locate  the  other  tip  so  that  the  crack  length 
could  be  measured.  This  method  was  quite  accurate  when  the  crack  was  stable  at  low 
load  levels.  But  as  the  fracture  load  was  approached,  the  tip  of  the  crack  whose  location 
was  first  determined  did  not  necessarily  remain  stable  while  the  microscope  was  moved 
to  locate  the  other  tip.  Furthermore,  in  some  cases,  the  tips  of  the  crack  were  not 
well  defined  and  precise  measurement  was  difficult.  However,  since  the  "critical” 
crack  length  for  each  specimen  was  determined  from  the  crack  extension  rate  which  was 
averaged  from  several  crack  measurements,  the  uncertainty  of  the  crack  measurements 
attributable  to  the  two  aforementioned  difficulties  was  considerably  reduced. 

3.  Constraint  introduced  by  the  side  grips  in  the  combined  tension  and 
shear  tests.  In  the  tests  which  followed  path  No.  3,  side  grips  were  required.  The 
constraint  of  the  side  grips  on  the  specimen  was  minimized  as  much  as  possible  by  the 
methods  outlined  in  sections  III  A2  and  III  B,  however  it  could  not  be  eliminated  during 
the  crack  extension  process. 

It  can  be  seen  that  in  the  combined  tension  and  shear  tests  which  followed  path 
No.  2  no  side  grips  were  required.  The  effect  of  the  constraint  of  the  grips  on  fracture 
strengths  of  the  specimens  can  be  examined  by  comparing  the  results  of  tests  which 
follow  path  No.  3  to  those  which  follow  path  No.  2.  In  Fig.  24,  the  data  which  are  pre¬ 
sented  as  open  points  were  obtained  from  tests  where  side  grips  were  used,  and  the  data 
presented  in  solid  points  were  obtained  from  tests  where  no  side  grip's  were  present. 
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In  the  middle  region,  where  both  types  of  experiments  were  performed  it  can  be  seen 
that  the  data  overlap  about  evenly.  Thus  it  can  be  concluded  that  the  constraint  of  the 
side  grips  was  properly  minimized,  and  its  effect  on  fracture  strength  cannot  be  de  - 
tected  from  the  experimental  results.  y 

4.  Buckling  of  the  specimen  under  shear.  Buckling  occurred  only  when 
shear  load  was  applied  to  the  specimen  with  side  grips  in  place.  When  buckling 
occurred,  the  fracture  strength  was  greatly  reduced,  but  when  shear  was  induced  by 
tests  which  followed  path  No.  2,  no  buckling  was  present.  Once  again,  the  data  in 
Fig.  24  indicated  that  there  was  no  significant  difference  between  the  data  obtained  by 
these  two  testing  methods,  and  it  can  be  concluded  that  buckling  was  successfully 
eliminated  by  the  constraint  system  employed  in  the  experiments*  when  shear  load  was 
applied  through  the  side  grips. 

5.  Absolute  similitude  was  not  maintained  between  the  theoretical  model 
in  the  analysis  and  the  model  used  for  experimental  investigation.  The  imposed 
conditions**  which  were  necessary  for  the  mathematical  analysis  of  the  crack  tip  stress 
distribution  will  be  compared  one  by  one  with  the  actual  conditions  of  the  experiments . 
First,  the  material  was  assumed  to  be  linearly  elastic.  The  specific  gravity  of  the 
specimens  varies  from  0. 126  to  0.220.  The  ultimate  tensile  strength  and  proportional 
limit  perpendicular  to  the  grain  for  balsa  wood  whose  average  specific  gravity  is  .  176 
ranges  from  118  to  170  psi  and  100  to  144  psi  respectively  ***.  No  data  on  the 
proportional  limit  for  pure  shear  was  available,  however  the  range  of  shearing 
strengths  is  from  298  to  360  psi.  In  the  experiments,  the  highest  nominal  stresses 
observed  were  112  psi  in  tension  and  258  psi  in  pure  shear.  Since  the  range  of  the 
strength  properties  is  relatively  large,  it  was  not  possible  to  examine  the  validity  of 
the  assumption  of  linear  elasticity  for  every  specimen,  but  it  is  clear  that  in  a  small 
region  around  the  crack  tips  the  stresses  exceeded  the  limits  of  elastic  behavior. 
Nevertheless,  from  the  range  of  the  values  given  above,  it  appeared  that  the  nominal 
stress  did  not  exceed  the  elastic  limit  with  the  exception  of  perhaps  a  few  specimens 
containing  very  short  cracks.  In  those  cases,  some  plasticity  correction  was  needed. 
However,  since  the  data  did  not  significantly  deviate  from  a  straight  line  of  slope  equal 
to  -1/2  in  the  log  u  vs  log  a  and  log  t  vs  log  a  diagrams  as  predicted  by  the  elastic 
analysis,  no  plasticity  correction  was  attempted. 

* 

**  This  is  described  in  section  III  A  2  and  shown  in  Fig.  8. 

***  These  conditions  are  stated  in  section  II  A. 

These  properties  are  taken  from  the  data  published  by  the  Balsa  Ecuador  Lumber 
Corporation,  500  F if th  Avenue ,  New  York  36,  New  York.  The  proportional  limits 
are  obtained  from  compression  tests. 
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Second,  the  plate  was  assumed  to  be  thin  and  to  extend  to  infinity  in  all  directions 
from  the  crack.  The  ratio  of  the  smallest  dimension  in  the  plane  of  the  plane  (the  width 

D  =4")  to  the  thickness  of  the  plate  t  was  -5-  =  64  which  indicated  that  plane  elasticity 
was  applicable.  Furthermore,  the  smallest  ratio  of  the  plate  width  W  (  at  the  line  of 

W 

the  crack)  to  the  crack  length  2a,  was  — jg  =3.6.  From  experience  with  isotropic 
plates,  it  was  known  that  edge  effects  of  the  plate  on  crack  extension  cannot  be  meas¬ 
ly 

ured  for  -j—  >  3.  It  is  therefore  reasonable  to  conclude  that  the  edge  effect  was  negli¬ 
gible  in  these  experiments . 

Third,  the  crack  was  assumed  to  be  straight  and  ’’sharp".  A  straight  crack 
was  cut  into  the  specimen  by  a  surgical  knife  blade  approximately  0.  02"  thick.  A  knife 
blade  whose  width  was  smaller  than  the  length  of  the  crack  was  inserted  in  the  center  and 
flexed  slightly  to  make  the  crack  tip  propagate  a  small  distance .  In  this  way  a  sharp 
natural  crack  was  introduced  in  the  specimen  which  corresponds  closely  to  the  theoreti¬ 
cal  assumption. 

Fourth,  the  crack  is  assumed  to  be  parallel  to  one  of  the  principle  planes  of 
elastic  symmetry.  In  the  experimental  program,  only  the  case  of  the  crack  parallel  to 
the  grain  of  the  wood  was  considered.  It  was  very  difficult  to  define  the  grain  of  the 
balsa  wood  precisely  and  consequently  some  error  was  introduced.  But  the  variations 
in  crack  orientation  were  within  +  2°  from  the  average  plane  of  elastic  symmetry  and  no 
provision  was  made  to  account  for  errors  thus  introduced. 

Fifth,  the  influence  of  environment  is  neglected.  If  the  environmental  conditions 
were  constant  for  all  of  the  experiments,  its  effects  could  have  been  neglected.  How¬ 
ever,  due  to  the  fact  that  there  was  no  means  of  controlling  the  environment  in  the 
laboratory  in  which  these  experiments  were  carried  out,  the  effects  of  temperature  and 
humidity  on  fracture  strength  were  not  accounted  for. 

Although  absolute  similitude  could  not  be  obtained  between  the  analytical  and 
experimental  model  without  unnecessarily  complicating  the  experiments  it  is  believed 
that  the  degree  of  similitude  achieved  was  sufficient  to  warrant  the  conclusions. 
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IV.  SUMMARY  AND  RECOMMENDATIONS  FOR  FURTHER  INVESTIGATION 


From  the  theoretical  analysis  of  the  stress  in  the  vicinity  of  the  tip  of  a  crack  in 

an  orthotropic  plate  it  was  found  that  the  stress  distribution  was  dependent  on  the  mate- 

-1/2 

rial  constants  and  that  stress  singularity  was  of  the  order  r  .  Also  the  theoretical 
analysis  indicated  that  stress  intensity  factors  =  oVa  and  k2  =  T-vCa  similar  to  those 
used  for  isotropic  materials  could  be  used  as  parameters  to  measure  the  strength  of  an 
orthotropic  plate  containing  a  crack.  Furthermore,  the  existence  of  a  functional  rela¬ 
tion  between  k^  and  k2  was  proposed. 

In  the  experimental  investigation  on  balsa  wood  plates,  the  stress  intensity  factor 
klc  and  k2c  were  found  to  be  constants  for  critical  combinations  of  crack  length  and 
loads  in  the  symmetrical  and  skew -symmetrical  cases  respectively.  Moreover,  the 
stress  intensity  factors  were  not  only  found  to  be  applicable  for  combined  symmetrical 
and  skew -symmetrical  loads,  but  they  were  related  by  the  function: 


k  1  k  ^ 

1  2 

(y2)  +  <  k1)  =1‘ 

lc  2c 


(60) 


From  the  experimental  results,  a  three  dimensional  representation  of  the  relationship 
among  the  variables  cr,  t  and  2a  is  given  in  Fig.  28.  In  this  figure,  the  surface  re¬ 
lating  a,  t  and  2a  can  be  considered  as  a  fracture  surface.  If  the  state  of  a  cracked 
orthotropic  plate  is  inside  the  tunnel  shaped  surface,  unstable  fracture  is  not  probable 
but  if  the  state  is  outside  of  this  surface  then  fracture  is  likely.  Thus,  it  can  be  seen 
that  fracture  mechanics  can  be  applied  to  orthotropic  materials. 

Since  the  analysis  deals  with  orthotropic  plates  of  only  one  layer,  the  analysis 
should  be  extended  to  the  more  realistic  case  of  multilayer  orthotropic  plates . 

The  experimental  program  was  designed  as  a  pilot  test  to  see  if  fracture  mechan¬ 
ics  could  be  applied  to  orthotropic  plates  as  suggested  by  the  theoretical  analysis  and  to 
explore  any  experimental  problems.  It  was  found  that  some  refinements  of  the  experi¬ 
mental  set  up  were  desirable.  A  tension-torsion  filament  wound  tube  specimen  appears 
attractive  for  future  work.  Experimental  investigation  of  multilayer  specimens  should 
be  considered  as  a  logical  extension  of  this  work  in  the  future  . 
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(a)  Crack  Tip  Shown  at  10X  (b)  Plate  Containing  Single  Crack  Subjected  to 

Symmetric  and  Skew-symmetric  Traction 
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Fig.  2  Superposition  of  Symmetric  Traction  on  Skew-symmetric  Traction 


(a)  Photograph  Illustrating  Irregularities 
of  Balsa  Wood  Specimen 


Fig.  4 


(b)Typical  Broken  Shear  Specimen 


Fig.  5  Testing  Machine 


(b)  Bottom  View 


6  Schematic  Diagram  of  Essential  Components 
of  Testing  Machine 


Fig.  8  Schematic  Diagram  of  Cross  Sectional  View  of  Experimental  Apparatus 
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Fig.  11  Typical  Crack  Extension  Patterns  of  Specimens  Subjected  to  Tension,  Path  No. 


320 


Fig.  12  Typical  Crack  Extension  Patterns  of  Specimen  Subjected  Fig.  13  Typical  Crack  Extension  Patterns  of  Speci 

to  Combined  Tension  and  Shear,  Path  No.  3  Subjected  to  Pure  Shear,  Path  No.  4 
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Critical  Half  Crack  Length  ,  a  ,  inches 
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Fig.  24  Interaction  Between  Critical  Symmetric  Stress  Intensity  Factor 
and  Critical  Skew-symmetric  Stress  Intensity  Factor  k2c 


Fig.  27  Effect  of  Variation  in  Specific  Gravity  of  Balsa  Wood 
on  Critical  Stress  Intensity  Factors 


Fig.  28  3~Dimensional  Representation  of  Experimenally  Determined  Fractured^  Surface 
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